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I  INTRODUCTION 

Stimulated  Raman  scattering  will  limit  the  maximum  laser  Intensity  that 
can  be  transmitted  through  the  atmosphere  If  the  divergence  of  the  beam  Is 
required  to  remain  constant.  To  establish  this  maximum  Intensity,  It  Is 
necessary  to  obtain  an  accurate  description  of  the  steady-state  Raman  gain 
coefficient  In  the  atmosphere.  This  report  describes  a  model  of  the  gain 
coefficient  for  the  Stokes  rotational  transitions  In  N2  under  atmospheric 
conditions.  This  model  Is  based  on  the  results  of  laboratory  measurements  for 
the  polarizability  anisotropy,  the  density  self-broadening  coefficients,  and 
the  O2  foreign-gas  density  broadening  coefficients.  These  measurements  are 
described  In  Appendices  A  and  B  and  are  used  with  a  published  atmospheric 
model^  to  calculate  Che  Raman  gain  coefficient  as  a  function  of  altitude  from 
0  to  100  km. 

It  should  be  emphasized  here  that  the  modeling  of  the  propagation  of  a 
high  Intensity  laser  pulse  through  a  Raman  medium,  such  as  the  atmosphere, 
requires  a  propagation  code  that  Includes  transient  effects,  off-axis  propaga¬ 
tion  for  the  generated  beams,  the  Inclusion  of  all  possible  Stokes  and 
anti-Stokes  orders  that  can  be  generated,  and  the  provision  for  a  variable 
pump  laser  bandwidth.  Codes  Including  some  of  these  effects  are  currently 
being  developed  in  several  laboratories.  The  model  of  the  steady-state  gain 
coefficient  presented  here  will  be  an  Important  input  to  these  codes. 


II  MODEL  DESCRIPTION 


A-  Atmospheric  Model 

In  the  atmospheric  model,  the  Raman  gain  calculation  uses  the  altitude- 
dependent  temperature,  N2  density,  and  O2  density  as  Input  data.  The 
atmospheric  model  used  here,  given  by  Banks  and  Kockarts,^  Is  for  Wallops  Island 
Virginia  (38°N).  Temperatures  and  densities  for  N2  and  O2  were  taken  from 
Table  3.1  of  Reference  1.  Densities  were  extrapolated  for  the  first  15  km 
(omitted  from  Table  3.1  of  Reference  1)  by  assuming  an  exponential  decay  with 
Increasing  altitude,  an  e~^  altitude  of  8  km,  and  relative  densities  of  78% 
and  21%  for  N2  and  O2,  respectively.  The  average  temperature  profiles  for 
summer  and  winter  are  shown  In  Figure  1.  Maximum  variations  are  Indicated  by 
horizontal  bars  for  summer  and  by  shading  for  winter.  Temperatures  for  the 
first  15  km  (also  omitted  from  Reference  1)  were  determined  from  the  average 
of  the  summer  and  winter  values  In  Figure  1. 


B.  Raman  Gain  Model 

The  steady-state,  plane  wave  Raman  gain  coefficient  Is  given  by"^ 
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where  f(5v)  Is  the  area-normalized  Raman  llneshape  and  6v  =  -  v^,  where 

Vp,g  Is  the  pump/Stokes  laser  frequencies,  Is  the  Raman  transition 
frequency,  and  AN  Is  the  population  difference  as  specified  In  Appendix  A.  At 
the  peak  (6v  «  0)  of  a  Lorentzlan-shaped  transition,  equation  (1)  from  above 
reduces  to  equation  (1)  of  Appendix  A.  This  study  uses  equation  (1)  to  calcu¬ 
late  the  Raman  gain  coefficient;  hence  the  spontaneous  scattering  cross 
section,  3a/3n,  and  the  llneshape,  f ( 6v) ,  must  first  be  determined. 
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TEMPERATURE  PROFILES  USED  DJ  THE  RAMAN  GAIN  CALCULATIONS, 
variations  are  indicated  by  horizontal  bars  for  summer  and 
shading  for  winter  (after  Ref.  1). 
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The  cross  section,  3a/3Q  Is  given  by  equation  (2)  of  Appendix  A.  This 
expression  contains  the  polarizability  anisotropy,  y,  which  is  given  by 
equation  (3)  of  Appendix  A  and  reproduced  here: 
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The  parameters  >  and  and  the  procedure  used  to  determine 

them  are  also  given  in  Appendix  A.  The  wavelength  dependence  of  equation  (2) 
is  shown  as  a  solid  line  in  Figure  2  along  with  several  experimental 
values.  Note  that  the  solid  line  is  not  a  direct  fit  to  the  data  points, 
but  is  derived  according  to  the  procedure  described  in  Appendix  A.  For  N2, 
the  J  dependence  of  the  polarizability  anisotropy  is  <  1%  for  0  <  J  <  20  and 
is  neglected.  Equation  (2)  gives  a  value  of  7.39  x  10~^^  cm^  for  Xg  •  350 
nm.  All  gain  results  presented  in  this  report  (except  those  in  Appendix  A)  are 


for  X, 


350  nm. 


The  accuracy  of  equation  (2)  can  be  checked  by  using  the  parameters 


v.  ,  a  ,  and  a  to  calculate  the  index  of  refraction,  n.  The  results 
XX  SR  SX  ^ 


can  then  be  compared  with  experimental  data  for  n^l: 
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This  comparison  is  made  in  Figure  3,  where  the  dashed  line  shows  our  current 
model  and  the  triangles  show  the  experimental  data  from  the  literature.  The 
agreement  is  excellent,  with  the  largest  discrepancy  for  (n  -  1)  being  about 
0.5%.  Figure  3  Illustrates  the  high  accuracy  reported  in  Appendix  A  for  the 
polarizability  anisotropy  of  equation  (2). 

The  population  difference,  AN,  was  calculated  using  a  N2  rotational  con¬ 
stant  of  1.99  cm  More  accurate  coefficients  can  be  found  in  Reference  8. 
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FIGURE  2  WAVELENGTH  DEPENDENCE  OF  THE  POLARIZABILITY  ANISOTROPY  DETER¬ 
MINED  FROM  LEAST  SQUARES  FITS  TO  THE  AVAILABLE  EXPERIMENTAL 
DATA  (INDIVIDUAL  POINTS).  See  Appendix  A  for  the  procedure 
used  to  generate  the  fit  (solid  line). 
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FIGURE  3 
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COMPARISON  OF  THE  EXPERIMENTAL  DATA  (A)  FROM  REF.  7  AND  THE 
RESULTS  OF  OUR  MODEL  (DASHED  LINE)  FOR  THE  INDEX  OF  REFRACTION,  n 


Pure  Voigt  lineshapes  were  used  for  all  results  presented  here.  The 
width  of  the  Lorentzian  contribution  to  the  Voigt  was  obtained  by  adding  the 
linewldths  due  to  self-broadening  and  foreign-gas  broadening  due  to  02*  The 
present  llnewldth  measurements  as  a  function  of  temperature  and  density  are 
summarized  in  Appendix  B.  Self-broadening  coefficients  not  given  in  Table  I 
of  Appendix  B  were  determined  by  assuming  a  linear  J  dependence.  Foreign-gas 
broadening  coefficients  for  temperatures  and  J’s  not  given  In  Table  V  of 
Appendix  B  were  determined  by  assuming  the  foreign-gas  values  are  85?!  of  the 
corresponding  self-broadening  values  In  Table  I  of  Appendix  B.  Tn  this  study, 
all  broadening  coefficients  from  Appendix  B  were  assumed  to  vary  linearly  with 
temperature. 

Only  Gaussian  laser  lineshapes  were  considered  for  this  model.  Thus,  the 
llnewldth  of  the  Gaussian  component  for  the  Voigt  profile  was  determined  from 
the  square  root  of  the  sum  of  the  squares  of  the  laser  and  Doppler  line- 
widths.  After  the  Lorentzian  and  Gaussian  linewldths  were  determined,  an 
area-normalized  Voigt  calculation  was  used  to  determine  f( 6v)  In  equation  (1). 


Ill  RESULTS 


A.  Altitude  Dependence 

Figure  4  shows  the  peak  Raman  gain  coefficient  as  a  function  of  altitude 

for  S(8).  For  comparison,  the  dotted  line  shows  the  gain  if  the  approxima- 
2  2  2 

tion  Av,p  =  AVq  +  is  used  with  equation  (1)  of  Appendix  A.  The  quantities 
5v.j.,  Avq,  and  Av^  total,  Gaussian,  and  Lorentzian  FWHM  linewidths, 

respectively.  The  two  quantities  AN  and  Av>p  that  produce  this  altitude 
profile  are  also  shown  as  functions  of  altitude  in  Figures  5  and  6.  The  Raman 
gain  of  Figure  4  Is  approximately  constant  below  40  km  and  drops  to  zero  above 
40  km.  This  occurs  because  the  peak  Raman  gain  depends  on  the  ratio  AN/Av.p. 
Below  40  km.  Figures  5  and  6  show  that  AN  and  Avq<  decrease  with  altitude  at 
about  the  same  rate,  yielding  a  constant  value  of  gain.  Above  40  km.  Figure  6 
shows  that  Av-j.  becomes  constant,  whereas  Figure  5  shows  AN  continues  to  de¬ 
crease;  thus,  the  Raman  gain  also  decreases.  The  Raman  llnewidth  is  constant 
at  high  altitudes  because  the  collislonal  broadening  becomes  negligible  com¬ 
pared  with  the  roughly  constant  Doppler  broadening. 

The  peak  Raman  gain  coefficients  for  the  strongest  Stokes  transitions  are 
shown  as  a  function  of  altitude  in  Figure  7.  The  number  indicated  next  to 
each  curve  is  the  lower  rotational  quantum  number  for  that  particular 
transition.  The  odd-numbered  J  transitions  are  only  half  as  strong  as  their 
adjacent  even-numbered  neighbors  and  are  omitted  from  this  Investigation.  The 
lowest-J  transitions,  S(0)  through  S(4),  have  gains  that  decrease  less  quickly 
with  altitude  and  are  shown  in  Figure  8,  where  the  vertical  scale  has  been 
magnified  by  a  factor  of  three  from  that  of  Figure  7,  The  integrated  gain  co¬ 
efficients  over  the  depth  of  the  atmosphere  (from  sea  level  to  100  km)  for  all 
these  transitions  are  tabulated  in  Table  1. 

Raman  gains  for  frequencies  other  than  those  at  the  peak  [5v  =  0  in 
equation  (1)]  of  the  Raman  transitions  were  also  calculated.  Figure  9  shows 
the  S(8)  gain  coefficient  as  a  function  of  altitude  for  detuning  values  of  0, 
10,  100,  and  800  MHz.  The  integrated  (over  altitude)  gain  coefficient  for 
S(8)  is  plotted  versus  detuning  in  Figure  10. 
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POPULATION  DIFFERENCE,  AN,  AS  A  FUNCTION  OF  ALTITUDE  FOR 
S(8)  IN  N2. 
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FIGURE  6  RAMAN  LINEWIDTH  (FWHM)  OF  VOIGT  LINESHAPE  AS  A  FUNCTION  OF 
ALTITUDE  FOR  S(8)  IN  N2 . 
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FIGURE  8 
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PEAK  RAMAN  GAINS  FOR  S(0)  THROUGH  S(4)  AS  A  FUNCTION  OF  AI.TITUDE. 
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S(8)  RAMAN  GAIN  AS  A  FUNCTION  OF  ALTITUDE  FOR  DETUNINGS  OF 
0,  10,  100,  AND  800  MHz  FROM  THE  LINE  CENTER. 


Seasonal  and  Laser  Llnewldth  Effects 


All  the  results  of  Figures  A  through  10  were  obtained  for  the  specific 
temperature  profile  of  Table  3.1  of  Reference  I.  The  effect  of  atmospheric 
temperature  fluctuations  is  summarized  in  the  four  curves  of  Figure  11.  The 
Raman  gain  of  S(8),  due  to  seasonal  variations  of  temperature,  is  shown  with 
two  gain  profiles  for  both  summer  and  winter.  The  curves  labeled  maximum  show 
the  gain  for  the  highest  temperatures  for  that  season,  whereas  those  labeled 
minimum  show  the  gain  for  the  lowest  temperatures.  Figure  11  shows  that 
maximum  variations  in  the  Raman  gain  due  to  temperature  changes  in  the 
atmosphere  will  be  about  10%  at  all  altitudes.  The  four  temperature  profiles 
used  in  the  calculations  of  Figure  11  are  taken  from  Figure  1. 

The  last  Raman  gain  parameter  discussed  here  is  the  laser  llnewldth.  All 
previous  results  in  this  paper  were  obtained  assuming  negligible  laser  line- 
widths.  Appreciable  Gaussian  laser  frequency  profiles  were  also  considered. 
The  effect  of  the  laser  bandwidth  on  the  gain  profiles  of  S(4)  and  S(6)  is 
shown  in  Figures  12  and  13,  respectively.  The  laser  linewidths  (in  llHz)  are 
indicated  next  to  each  curve.  Integrated  gains  for  all  the  cases  in  both  of 
these  figures  are  listed  in  Table  2.  Figures  12  and  13  show  that  in  the 
absence  of  phase  locking  between  the  pump  and  Stokes  waves, large  laser 
bandwldths  can  effectively  reduce  the  gain  length. 

C.  Voigt  Approximation 

In  this  study  the  Raman  llneshape  was  always  modeled  as  a  Voigt  func¬ 
tion.  In  this  section,  a  simple  approximation  (requiring  less  computation 
time  than  an  exact  Voigt  calculation)  is  presented  as  an  alternative  to  calcu¬ 
lating  the  Voigt  function.  This  approximation  is  valid  when  calculating  the 
peak  Raman  gain  coefficient,  but  is  not  applicable  for  frequencies  off  the 
line  center. 

The  procedure  uses  equation  (1)  to  calculate  the  peak  Raman  gain.  The 
llneshape,  f(6v),  is  approximated  as 
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FIGURE  11  S(8)  PEAK  RAMAN  GAIN  VARIATIONS  FOR  MAXIMUM  VARIATIONS  IN 

ATMOSPHERIC  TEMPERATURES . 
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FIGURE  12  S(4)  GAIN  AS  A  FUNCTION  OF  ALTITUDE  FOR  LASER  LINEWIDTHS  OF 

0,  50,  250,  AND  500  MHz.  This  calculation  neglects  the  phase¬ 
locking  of  the  pump  and  Stokes  waves  for  collinear  propagation.  If 
the  Stokes  wave  builds  up  from  noise  (not  seeded),  the  steady-state 
Raman  gain  will  be  that  of  the  curve  for  0  MHz  (assuming  dispersion 
is  negligible). 
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FIGURE  13  S(6)  GAIN  AS  A  FUNCTION  OF  ALTITUDE  FOR  LASER  LINEWIDTHS  OF 

0,  50,  250,  AND  500  MHz.  This  calculation  neglects  the  phase- 
locking  of  the  pump  and  Stokes  waves  for  collinear  propagation.  If 
the  Stokes  wave  builds  up  from  noise  (not  seeded),  the  steady-state 
Raman  gain  will  be  that  of  the  curve  for  0  MHz  (assuming  dispersion 
is  negligible) . 


where  Avq  and  are  the  total  Gaussian  and  Lorentzian  FWHM  linewldths.  The 
difference  between  the  approximation  of  equation  (4)  and  an  exact  Voigt 
calculation  Is  Illustrated  In  Figure  14,  where  the  Raman  gain  for  both  cases 
Is  plotted  as  a  function  of  the  ratio  of  the  Lorentzian  and  Gaussian  contribu¬ 
tions  to  the  llnewldth.  The  fractional  difference  of  f(6v=0)  for  the  two 
methods  Is  also  shown,  and  has  a  maximum  value  of  about  3%.  Results  In  Figure 
14  are  obtained  with  a  constant  Gaussian  FWHM  of  5  MHz,  corresponding  to  the 
Doppler  broadening  of  S(8)  at  room  temperature. 


LORENTZIAN  WIDTH/GAUSSIAN  WIDTH 


JA-7193-71 


FIGURE  14  COMPARISON  OF  f(6v=0)  USING  THE  EXACT  VOIGT  LINESHAPE  (DASHED 
CURVE)  AND  THE  APPROXIMATION  GIVEN  BY  EQUATION  (4)  (DOTTED 
CURVE).  The  fractional  difference  (solid  curve)  is  also  shown. 
The  variable  along  the  horizontal  axis  is  the  ratio  of  the 
Lorentzian  llnewidth  to  the  Gaussian  linewidth  and  the  Gaussian 
FWHM  is  held  constant  at  5  MHz. 


IV  SUMMARY 


The  Raman  gain  of  the  Stokes  rotational  lines  In  N2  has  been  modeled  as  a 

function  of  altitude  for  the  lower  100  km  of  the  atmosphere.  Calculations 

were  completed  for  the  even-J  transitions,  S(0)  through  S(16).  The  strongest 

lines,  S(6)  through  S(10),  have  gain  coefficients  that  are  approximately  con- 
2  —1  —1 

stant  at  1  cm  MW  km  over  the  lower  40  km  and  then  decrease  to  zero  over 
the  40-80  km  region.  Effects  on  the  Raman  gain  are  described  for  seasonal 
temperature  fluctuations  In  the  atmosphere,  laser  llnewldths,  and  detuning 
from  the  line  center  of  the  Raman  transition. 

The  uncertainties  In  the  Raman  gain  coefficients  reported  here  depend  on 
how  well  the  temperature  of  the  gain  medium  is  known.  If  the  temperature  Is 
accurately  known  (e.g.,  laboratory  conditions),  the  uncertainties  In  the  gain 
coefficients  are  limited  by  the  uncertainties  in  the  llnewldths  and  the 
polarizability  anisotropy.  Both  of  these  uncertainties  are  about  1-2%;  thus, 
the  gain  coefficient  uncertainties  are  about  5%.  In  the  atmosphere,  where  the 
temperature  Is  not  well  known,  the  uncertainties  In  the  gain  coefficients  are 
dominated  by  temperature  uncertainties.  Therefore,  gain  coefficient  uncer¬ 
tainties  are  about  10%,  as  Illustrated  In  Figure  11. 
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Appendix  A 


TEMPERATURE  AND  WAVELENGTH  DEPENDENCE  OF  THE  ROTATIONAL 
RAMAN  GAIN  COEFFICIENT  IN  N2 

G.  C.  Herring,  Mark  J.  Dyer,  and  William  K.  Bischel 
Chemical  Physics  Laboratory 
SRI  International,  Menlo  Park,  CA  9402 5 

ABSTRACT 

The  temperature,  wavelength  and  J  dependence  of  the  rotational  Raman  gain 
coefficient  has  been  determined  for  the  S  branch  transitions  in  H2*  First, 
the  temperature  dependence  (80-300K)  of  the  density-broadening  coefficients 
was  measured  using  stimulated  Raman  spectroscopy.  Second,  the  wavelength 
dependence  (250-600nm)  of  the  polarizability  anisotropy  was  emperically 
determined  by  fitting  to  the  most  recent  experimental  data.  These  two  results 
were  then  used  to  calculate  the  rotational  Raman  gain  coefficients  with 
accuracies  estimated  at  5%.  At  room  temperature,  the  high  density  limit  of 
the  steady-state  gain  coefficient  of  the  strongest  line,  S(10),  is  4.8  *  10“^^  cra/W 
for  a  Stokes  wavelength  of  568  nm. 
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There  are  many  applications,  such  as  laser  fusion,  that  require  the  pro¬ 
pagation  of  high  Intensity  lasers  through  the  atmosphere.  Recently,  there  has 
been  renewed  Interest  In  the  nonlinear  optical  processes  that  limit  the  maxi¬ 
mum  transmitted  laser  Intensity  for  a  given  propagation  distance.  Several 
different  nonlinear  effects  may  be  Important  and  some  of  these  were  summarized 
by  Zuev^  in  1982.  However,  It  has  only  recently  been  suggested  that  rotatio¬ 
nal  stimulated  Raman  scattering  (RSRS)  in  N2  may  ultimately  limit  the  maximum 
Intensity  that  can  be  propagated  through  air  If  the  frequency  and  divergence 
of  the  laser  are  to  remain  unchanged.  This  suggestion  Is  supported  by  recent 
experiments*  at  Lawrence  Livermore  National  Laboratory,  where  RSRS  In  N2  wns 
observed  when  the  NOVA  fusion  laser  was  propagated  through  100  m  of  air  at 
Intensities  In  the  range  of  a  few  GW/cm^. 

This  observation  has  reinforced  the  need  for  accurate  modeling  of  stlmu- 
lated  Raman  scattering  In  air.  Unfortunately,  the  existing  uncertainties**^ 

In  the  steady-state  gain  coefficient  in  N2  limit  the  usefulness  of  any 
potential  model.  In  particular,  no  data  exist  for  the  temperature  dependence 
of  the  gain  coefficient,  an  Important  consideration  in  any  model.  In  this 
letter  we  present  new  experimental  results  for  the  temperature  dependence  of 
the  Raman  gain  coefficient  for  the  most  Important  S  branch  rotational  transi¬ 
tions  In  N2. 

Direct  measurement  of  the  gain  coefficient  requires  an  accurate  measure¬ 
ment  of  the  laser  Intensity  In  the  interaction  region.  For  low  gain  Raman 
systems  such  as  N2,  this  Is  extremely  difficult  since  a  focused  geometry  Is 
usually  required  to  obtain  adequate  signal  strengths.  Alternatively,  It  has 
been  recently  experimentally  demonstrated^  that  the  gain  coefficient  can  be 
accurately  calculated  If  the  llneshape  function  for  the  Raman  transition  and 
the  polarlzablillty  anisotropy  are  known  precisely.  We  have  chosen  this 


alternative  method  for  determining  the  rotational  Raman  gain  coefficients  in 


N2«  Using  stimulated  Raman  gain  spectroscopy,^  we  have  measured  the  density 
and  temperature  dependence  of  the  Raman  llnewidths  for  the  Stokes  branch  rota¬ 
tional  transitions  in  N2>  In  addition,  the  relative  scattering  cross  sections 
were  measured  for  these  same  rotational  lines.  Using  the  results  of  these 
measurements  and  previously  published  values”  of  the  polarizability 
anisotropy,  we  have  calculated  the  high-d^nsity  Raman  gain  in  the  temperature 
region  80-300  K  for  four  of  these  rotational  transitions. 

For  a  Lorentzian  llneshape  function,  the  peak  plane  wave,  steady-state 
Raman  gain  coefficient  is^^ 
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where  Vg  is  the  Stokes  frequency  in  hertz,  Xg  is  the  Stokes  wavelength  in  the 
medium,  Av  is  the  llnewldth  (FWHM)  in  hertz,  and  AN  is  the  population 
difference,  which  is  equal  to  N( J)-[N( J' )(2J+1) /(2J'+1 ) ) .  N(J)  is  the 
Boltzmann  statistical  population  density  in  a  state  with  rotational  quantun 
number  J.  The  rotational  Raman  scattering  cross  section  for  incident  and 
scattered  waves  linearly  polarized  in  the  same  direction  is  given  hy^^ 

n 

da  2  ^  \  \  (J+1)  (J+2)  2 
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where  y  is  the  polarizability  anisotropy.  The  llnewldth  in  Eq.  (1)  varies 
with  J,  density  and  temperature,  while  y  varies  with  pump  laser  frequency. 
Thus,  it  is  necessary  to  know  this  parametric  dependence  in  order  to 
accurately  calculate  the  Raman  gain. 


The  experimental  setup  Is  a  quasl-cw  stimulated  Raman  spectrometer 
similar  to  that  In  Ref.  5.  The  cw  probe  laser,  which  provides  150  mW  of  power 
In  the  Interaction  region.  Is  a  single-mode  Kr"*”  Ion  laser  operating  at 
568  nm.  The  pump  laser  Is  a  coherent  699-29  single-mode  cw  ring  dye  laser 
that  Is  pulse  amplified  using  a  Quanta-Ray  PDA-1 .  At  565  nm,  the  tunable 
output  of  the  pump  laser  consists  of  3-mJ,  10-nsec  (FWHM)  pulses  at  10  Hz. 

The  frequency  resolution  of  this  system,  limited  by  the  llnewldth  of  the 
pulsed  laser.  Is  approximately  100  MHz.  A  250-ijm-dlameter  pinhole  Is  used  to 
Improve  the  spatial  mode  of  the  pump  beam  and  to  verify  that  the  beam  does  not 
move  when  the  dye  laser  Is  tuned  over  wavelength  Intervals  of  100  A.  The  pump 
and  probe  beams  were  focused  and  crossed  (angle  of  1  degree)  at  the  center  of 
a  gas  sample  cell  using  a  40-cm  focal  length  lens.  After  the  cell,  stray  pump 
light  was  Isolated  from  the  probe  beam  with  a  second  pinhole  while  an  AC-coupled 
photodiode  was  used  to  monitor  the  probe  Intensity.  The  Raman  gain  signal, 
which  appears  as  a  10-nsec  pulse  on  the  probe  laser  Intensity,  Is  averaged 
with  a  gated  Integrator  and  stored  on  disk  for  later  analysis. 

Density-broadening  coefficients  were  determined  for  the  S(6),  S(8), 

S(10),  and  S(12)  rotational  transitions  In  N2.  The  data  were  taken  at  three 
temperatures  (295  K,  195  K,  and  80  K)  for  densities  spanning  the  range 
0.01-2  amagat.  The  results  of  these  measurements  are  shown  In  Figure  1  and 
listed  In  Table  I.  The  values  for  the  density-broadening  coefficient  at  room 
temperature  are  In  agreement  with  ab  Initio  theory^^  (solid  line)  and  are  In 
good  agreement  with  recent  llnewldth  measurements  using  spontaneous 
scattering.  Simultaneous  recording  of  the  Raman  signals  at  two  different 
pressures  also  allowed  the  observation  of  density  shifts  of  the  Raman  transi¬ 
tion  frequency.  We  have  determined  that  the  density  shifts  are  less  than 
150  MHz/amagat  for  all  rotational  lines  and  temperatures  studied,  "^hese 


broadening  and  shift  measurements  will  be  described  In  greater  detail  In  a 
future  publication.^^ 

The  ground  state  polarizability  anisotropy  Is  defined  as 
Y  “  Oj  -  where  (i,l)  Indicate  the  perpendicular  and  parallel  components, 
respectively,  of  the  polarizability  a.  We  can  model  the  wavelength  dependence 
of  Y  using  an  empirical  formula® 


(3) 


where  v  Is  the  pump  laser  frequency,  Og .  ,  Is  the  static  or  DC  polarizability, 

and  I  la  an  effective  Intermediate  state.  The  ^  were  determined  first 

by  fitting  Eq.  (3)  to  the  results  of  recent  ab  Initio  calculations^  for  and 

The  Ogj  and  Og^  were  then  obtained  by  fitting  Eq.  (3),  with  I  fixed, 

to  recent  experimental  determinations® of  a.  and  a,  to  obtain  a_ .  and 

II  s  j 

Ogj^.  The  N2  parameters  derived  using  this  procedure  are: 

-  1.260  X  10^  cm"^,  -  1.323  x  10^  cm”^,  ^  -  2.200  x  lO"^^  cm^, 

and  a  ,  ■  1.507  x  10  cm^. 
si 

In  Fig.  2,  the  wavelength  dependence  of  y( v) ,  from  Eq.  (3),  Is  shown 
along  with  several  experimental  values.  We  have  used  data  from  only  Kefs.  6, 

8  and  9  to  determine  the  Og^^  I  because  these  measurements  have  substantially 
smaller  uncertainties  than  the  other  Refs,  cited  in  Fig.  2.  Eqs.  2  and  3 
yield  cross  sections  that  are  In  agreement  with  direct  measurements.^®*^® 

In  N2,  the  J  dependence  of  the  polarizability  anisotropy,  Yi  is 
negllglble^^ for  J  ^  20.  We  have  verified  this  conclusion  by  determining 
the  relative  cross  sections  from  our  stimulated  Raman  llneshape  measure¬ 
ments.  Thus  we  have  used  Eq.  (3)  to  determine  the  polarizability  anisotropy 


for  all  transitions  studied  In  this  letter 


The  reliability  of  Eq.  (3)  can  be  checked  by  using  the  parameters  a 


si,  n 


and  I  to  calculate  the  Index  of  refraction,  n;  the  result  can  then  he 
compared  with  experimental  data.  The  expression  for  n-1  Is 


o  2  , 

2  ,  ,  ^  n-1 

3  (n  -  1)  - 

n  +2 


4xN 

9 


/A 

‘sil  2  2r  “si  2 


Ai"  '' 


AT  " 


(4) 


where  N  is  the  density.  Over  the  wavelength  range  250-600  nm,  Eq.  (4)  agrees 

1  q 

with  experimental  data  to  better  than  0.5%,  verifying  a  high  accuracy  for 
Y(v)  of  Eq.  (3). 

We  have  used  Eqs.  1-3  and  our  linewldth  data  to  calculate  the  Raman  gain 
coefficient  for  the  rotational  transitions  as  a  function  of  temperature  In  the 
limit  where  density  broadening  is  the  dominant  contributor  to  the  linewldth. 
For  N2,  this  hlgh-denslty  limit  corresponds  to  densities  greater  than 
0.01  amagat.  The  results  are  summarized  In  Table  I,  while  the  temperature 
dependence  Is  Illustrated  In  Figure  3.  These  results  were  obtained  bv  assum¬ 
ing  the  rotational  llnewldths  are  linear  with  temperature  and  the  polariz¬ 
ability  anisotropy  Is  constant  with  J. 

The  accuracy  of  the  gain  coefficients  reported  here  Is  limited  by  the 
uncertainties  In  the  llnewldths  and  the  polarizability  anisotropy.  Our  line- 

width  measurements  have  uncertainties  of  2%  (one  standard  deviation)  as  given 

14 

by  the  linear  fits  used  to  determine  the  broadening  coefficients.  The 
uncertainty  In  the  polarizability  anisotropy  was  estimated  to  be  1.5%  from  the 
differences  In  the  results  of  different  Investigators.  Thus  the  accuracy  of 
the  gain  coefficients  reported  here  Is  approximately  5%. 

We  can  compare  the  pure  N2  gain  coefficients  reported  here  with  the  air 
coefficients  reported  by  Heneslan  et  al.,^  since  we  have  also  measured^^  the 
forlegn  gas  broadening  coefficients  of  O2.  After  correcting  for  N2  and 


densities,  Xg,  and  the  wavelength  dependence  of  the  polarizability  anlstropy 
(Eq.  3),  we  obtain  a  N2  gain  coefficient  of  2.0  x  10“^^  cm/W  for  one 
atmosphere  (80%/20%  mixture  of  N2/®2^  ®*-  ^s  “  This  Is  25%  smaller 

than  the  value  of  Heneslan  et  al.^,  but  conslstant  with  their  error  bar  of  a 
few  tens  of  per  cent. 

In  conclusion,  we  have  made  llnewldth  and  cross  section  measurements  that 
will  facilitate  the  modeling  of  stimulated  rotational  Raman  scattering  In 
N2.  Llnewldths  are  linear  with  density  for  densities  from  0.01  to  2.0  amagat 
and  are  also  linear  with  temperature  for  temperatures  from  80  K  to  300  K. 
Within  the  accuracy  (±  5%)  of  our  relative  measurements,  the  polarizability 
anisotropy  was  found  to  be  Independent  of  J  for  S(2)  through  S(16).  The  new 
llnewldth  data  presented  here  can  now  be  used  to  accurately  (±5%)  predict  the 
rotational  Raman  gain  over  a  wide  range  of  temperatures  and  densities  In  N2« 

In  addition  to  rotational  Raman  scattering,  vibrational  Raman  scattering 
Is  appreciable  for  the  Intensities  used  In  this  work.  A.  room  temperature 
study  has  already  been  reported,  and  future  work  in  this  area  should  include 
an  Investigation  of  the  temperature  dependence  of  the  gain  for  the  vibrational 
transition  In  N2. 

This  work  was  supported  by  the  Defense  Advanced  Research  Projects  Agency 
under  Contract  N00014-84-C-0256,  through  the  Office  of  Naval  Research. 
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Figure  2. 


Figure  3. 


Measured  temperature  dependence  of  the  density-broadening 
coefficient  for  Stokes  branch  rotational  Raman  transitions  in 
N2.  The  solid  line  indicates  the  most  recent  room  temperature  ab 
initio  calculations  given  in  Ref.  12. 

Wavelength  dependence  of  the  1^2  polarizability  anisotropy.  Solid 
line  is  given  by  Eq.  (3),  which  was  derived  with  a  fit  to  the  data 
of  Refs.  6,  8,  and  *).  Data  of  Refs.  2  and  IS  were  not  used  in  the 
fit,  but  are  shown  for  comparison. 

Temperature  dependence  of  the  plane  wave,  steady-state  Raman  gain 
in  N2  calculated  from  Eq.  1  for  the  hlgh-denslty  limit  (p  >  .01 
amagat).  The  pump  and  Stokes  polarizations  are  linear  and 


parallel. 


Table  I.  Temperature  dependence  of  rotational  Raman  llnewldths,  fractional 
population  differences,  and  gain  coefficients  for  S  branch  in 
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ABSTRACT 

The  temperature  and  density  dependence  of  the  rotational  Raman  llnevldths 
and  llneshlfts  for  the  diatomic  molecules  N.^  and  H2  has  been  measured  using 
stimulated  Raman  gain  spectroscopy.  Room  temperature  results  for  the  density 
broadening  coefficients,  B,  In  N2  are  compared  with  ab  Initio  calculations, 
while  the  low  temperature  (80  K  and  195  K)  results  are  the  first  to  be 
reported.  Values  of  y  In  the  relation,  B  «  T^,  were  determined  from  fits  to 
the  data  for  self-broadening  in  N2  yielding  0.25  <  y  <  0.39.  Foreign-gas  (O2) 
broadening  coefficients  for  the  N2  transitions  were  also  measured  and  found  to 
be  10-15Z  smaller  than  the  self-broadening  coefficients.  Our  llnewldth  and 
lineshlft  results  for  H2  are  in  general  agreement  with  previous  measurements. 
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I .  INTRODUCTION 


Observation^  of  stimulated  rotational  Raman  scattering  in  air  with  the 
NOVA  laser  at  Lawerence  Livermore  National  Laboratory  has  generated  interest 
in  the  modeling  of  the  Raman  gain  in  N2«  Accurate  knowledge  of  the  Raman 
llnewldths  is  necessary  for  any  useful  model  because  the  peak  Raman  gain  is 
Inversely  proportional  to  the  llnewldth.  Results  of  previous  studies^”^  of 
the  room  temperature  rotational  Raman  llnewldths  in  N2  typically  differ  by 
30%,  while  results  at  other  temperatures  have  not  been  reported.  The  main 
purpose  of  this  work  is  to  measure  the  temperature  dependence  of  the  line- 
widths  and  llneshlfts  for  the  rotational  Raman  transitions  (S-branch)  in  N2> 

This  work  differs  from  previous  measurements  of  rotational  Raman  line- 
widths^"^  in  three  aspects.  First,  ve  have  measured  the  temperature  depen¬ 
dence  of  the  line  broadening  over  the  range  80-300  K.  Previous  work  in  N2  has 
been  only  for  room  temperature.  Second,  the  densities,  used  in  this  work  (0.01- 
2.0  amagat)  are  two  orders  of  magnitude  lower  than  those  (1.0-100  amagat)  used 
in  previous  measurements.  Complications  due  to  overlapping  of  adjacent  lines 
are  eliminated  by  working  at  lower  densities  where  the  Raman  lines  are  typi¬ 
cally  separated  by  100  llnewldths.  Finally,  all  previous  measurements  have 
used  spontaneous  Raman  scattering,  whereas  we  have  used  stimulated  Raman 
scattering,  which  has  a  frequency  resolution  several  orders  of  magnitude 
higher  than  its  spontaneous  counterpart. 

This  paper  reports  experimental  results  for  line  broadening  and  llneshlft 
coefficients  for  Stokes  rotational  Raman  transitions  in  N2  and  H2*  Results 
are  presented  for  the  even  J  transitions,  S(2)  through  S(16),  at  temperatures 
of  80  K,  195  K,  and  295  K  in  N2  and  for  S(0)  and  S(l)  at  temperatures  of  80  K 
and  295  K  in  H2.  The  foreign  gas  (O2)  broadening  and  shifts  of  the  four 


strongest  N2  transitions  have  also  been  measured  at  195  K  and  295  K.  The  room 
temperature  broadening  coefficients  In  N2  are  compared  with  ab  Initio  calcula¬ 
tions,  «rtiereas  the  low  temperature  N2  self-broadening  and  the  O2  foreign  gas 
results  are  the  first  to  be  reported.  The  present  temperature  dependence  of 
N2  is  compared  to  that  found  In  studies  of  CO  and  C02*  Finally  the  present  H2 
results  are  compared  with  previous  measurements. 

II.  EXPERIMENTAL 

A.  Apparatus 

An  overview  of  the  experimental  setup  Is  shown  In  Fig.  1.  This  quasl-cw 
stimulated  Raman  spectrometer  Is  similar  to  that  developed  by  Esherick  and 
Owyoung  and  consists  of  a  cw  probe  laser;  a  tunable,  pulsed  pump  laser;  a 
pair  of  gas  (sample  and  reference)  cells;  and  fast  photodiodes  to  detect  the 
Induced  Raman  gain  on  the  probe  beams. 

The  probe  laser  is  a  Kr'*'-lon  laser,  operating  at  568  nm  and  forced  into 
single-mode  operation  with  the  Insertion  of  a  temperature-stabilized  etalon 
Into  the  cavity.  This  laser  has  a  time-averaged  llnewldth  less  than  30  MHz. 
The  probe  beam  was  split  Into  two  equal  intensity  beams  so  that  Raman  signals 
could  be  simultaneously  obtained  from  two  separate  gas  samples.  After  passing 
through  their  respective  sample  cells,  both  of  the  cw  probe  beams  were  chopped 
at  10  Hz  with  a  mechanical  chopper  (not  shown  in  Fig.  1).  This  resulted  In 
170  ps  square-wave  pulses  that  were  incident  on  the  photodiodes.  Chopping  of 
the  probe  increases  the  saturation  intensity  for  the  photodiodes.  The  average 
power  of  the  probe  beams  was  150  mW  in  both  the  sample  and  reference  cells. 


The  pump  laser  was  a  tunable,  single-mode  (1-MHz  llnevidth)  cw  dye  laser 
that  was  pulse  amplified  using  a  Quanta-Ray  PDA-1.  The  cw  dye  laser  Input  to 
the  PDA-1  was  typically  400  mW.  The  PDA-1  was  pumped  with  the  532-nm  output 
(130  mJ  per  pulse)  from  a  frequency  doubled  TAG  laser.  At  565  nm,  the  tunable 
output  of  the  PDA-1  consisted  of  2-mJ,  lO-nsec  (FWHM)  pulses  at  10  Hz.  The 
llnewldth  of  the  pump  laser  was  100  MHz,  about  a  factor  of  two  larger  than  the 
Fourier  transform  limit  for  a  10-ns  pulse.  Peak  pump  Intensities  at  the  focal 
planes  were  estimated  to  be  5  GH/cm^  In  the  sample  cell  and  70Z  of  that  In  the 
reference  cell. 

The  pump  and  probe  beams  were  mode-matched  and  then  were  focused  and 
crossed  In  the  gas  cells  using  40-cm  focal  length  lenses.  The  crossing  angle 
was  1.0  degree  In  both  of  the  cells.  A  series  of  dlaphrams  placed  around  the 
probe  beams  was  used  to  eliminate  scattered  pump  light.  Both  probe  beams  were 
monitored  with  reversed  biased  photodiodes  that  had  rise  times  of  less  than 
1.0  ns,  and  were  terminated  with  50-Q  resistors.  The  total  probe  power  Inci¬ 
dent  on  the  photodiodes  during  the  170- ps  probe  laser  pulses  was  10  mW.  Capa¬ 
citors  (300  pf)  were  used  to  block  the  170- ps  pulses  and  to  pass  the  10-ns 
Raman  signals.  These  signals  were  amplified  by  a  factor  of  100  and  then 
averaged  with  gated  Integrators.  Output  time  constants  for  the  Integrators 
and  total  scan  times  were  usually  0.5  s  and  4  min,  respectively.  Occasio¬ 
nally,  for  weak  transitions,  time  constants  and  scan  times  were  Increased  to 
2.5  s  and  10  min,  respectively.  After  the  pump  beam  exited  the  reference 
cell.  It  was  also  monitored  with  a  third  photodiode.  This  signal  was  used  to 
normalize  the  two  Raman  signals  for  slow  pump  power  variations  during  the 
scans. 

The  sample  cell  was  constructed  by  attaching  a  pair  of  3.0-cm-dlameter 
brewster  window  extensions  to  a  drllled-out  10  cm^  cube  of  aluminum.  One  face 


of  the  aluminum  cube  is  in  contact  vlth  a  second  reservoir  that  can  be  filled 


with  either  liquid  N2  or  an  acetone  dry-ice  bath.  The  cooled  sample  cell  was 
insulated  with  a  vacuum,  and  the  temperature  was  monitored  with  a  thermocouple 
that  was  Imbedded  In  the  aluminum  cube.  The  reference  gas  cell  was  kept  at 
room  temperature  and  a  pressure  of  10  torr  for  N2  scans  and  200  torr  for  most 
of  the  H2  scans. 

The  temperature  dependence  of  the  linewldths  and  llneshlfts  was  measured 
by  taking  measurements  at  three  temperatures:  80  K,  195  K,  and  295  K.  At  each 
temperature  the  llneshape  was  recorded  for  several  densities  over  the  range 
0.01-2.0  amagat.  For  the  low  temperature  runs,  there  was  a  time  period  rang¬ 
ing  from  a  few  to  15  minutes  between  scans  with  different  densities,  so  that 
any  gas  added  between  scans  always  had  sufficient  time  to  come  to  thermal 
equilibrium  with  the  cell.  For  each  set  of  runs  at  constant  temperature,  the 
density  was  randomly  changed  rather  than  constantly  adding  or  eliminating  gas 
from  the  sample  cell  for  each  succeeding  run. 

The  pressure  In  the  experimental  cells  was  measured  using  a  10,000-torr 
MKS  Baratron  gauge,  and  the  temperature  was  measured  using  a  Chrome 1-Alumel 
thermocouple  gauge  referenced  to  273  K.  The  density  for  both  N2  and  H2  was 
calculated  using  the  perfect  gas  law  since  the  density  calculated  using  the 

Q 

second  vlrlal  coefficient  deviates  from  this  calculation  by  less  than  1%  at 
the  densities  and  temperatures  used  In  these  experiments. 

B.  Data  Analysis 

Since  the  laser  llnewldth  was  comparable  to  the  low  density  Raman 
llnewldth.  It  was  Important  to  account  for  this  Instrumental  width  in  the  data 
analysis.  Although  the  pump  laser  frequency  profile  was  not  directly 
measured.  It  was  Indirectly  determined  by  scanning  the  laser  over  a  narrow 


Raman  line  (e.g.,  S(0)  In  H2  at  5  tort).  This  transition  has  negligible 
(1  MHz)  colllslonal  broadening  compared  to  100  MHz  of  Doppler  broadening.  The 
profile  for  this  transition  was  found  to  fit  extremely  well  to  a  pure 
Gaussian.  Since  the  convolution  of  two  Gausslans  results  In  another  Gaussian, 
we  concluded  that  the  pump  laser  profile  could  be  approximated  as  a 
Gaussian.  Fitting  Voigt  profiles  to  low  density,  narrow  N2  llneproflles  with 
known  Doppler  and  Lorentzlan  contributions  confirmed  this  conclusion.  The 
probe  laser  llnewldth  was  assumed  to  be  negligible  when  convolved  with  the 
pump  llnewldth;  thus  all  the  data  presented  here  were  analyzed  assuming  a 
Gaussian  profile  for  the  Instrumental  width. 

The  fitting  procedure  used  to  determine  the  Lorentzlan  portion  of  the 
llneshapes  was  different  for  the  N2  and  H2  data.  We  first  describe  the  fitt¬ 
ing  procedure  for  the  N2  transition.  For  high  densities  in  N2,  Doppler 
broadening  (10  MHz)  and  the  laser  llnewldth  (100  MHz)  can  be  neglected.  Thus 
sample  cell  llneshapes  for  high  densities  were  fit  to  pure  Lorentzians  to  get 
a  first  approximation  to  the  Lorentzlan  width.  The  low  density  reference  cell 
data  were  then  fit  to  Voigt  llneshapes  with  a  fixed  Lorentzlan  (determined 
from  the  high  density  fits)  and  variable  Gaussian  components.  This  Gaussian 
component,  representing  the  pump  laser  profile  for  that  particular  scan,  was 
then  used  In  a  final  Voigt  fit  to  the  sample  cell  profile.  Any  fluctuations 
In  the  pump  laser  llnewldth  from  scan  to  scan  were  accounted  for  with  this 
procedure.  The  change  In  Lorentzlan  llnewldth  determined  from  the  first 
Lorentzlan  fit  and  the  final  Voigt  fit  was  typically  a  few  percent  and  further 
Iterations  were  not  necessary. 

For  the  H2  transitions  studied,  experimental  constraints  did  not  allow 
the  density  In  the  cell  to  be  larger  than  3  amagat.  Therefore  the  approximate 
density  broadening  coefficient  could  not  be  determined  from  the  high  density 
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data.  Hence,  the  above  procedure  was  not  applicable.  Instead,  all  of  the  H2 
llneshapes  were  fit  to  Voigt  profiles  with  freely  adjustable  Lorentzlan  and 
Gaussian  components. 


III.  RESULTS 


A.  Line  Broadening  In  N2 

An  example  of  the  S(10)  Raman  profile  In  N2  at  195  K  Is  shown  In  Fig  2. 

The  narrower  profile  was  taken  at  10  torr  and  the  wider  profile  was  taken  at 

400  torr.  Individual  dots  represent  the  data,  and  the  sol  d  line  represents 
the  fit  for  a  Voigt  llneshape,  which  Is  dominated  by  Its  Lorentzlan  component 
for  the  400-torr  data.  A  slight  asymmetry  Is  observable  In  the  data  of  Fig.  2 

and  was  characteristic  of  most  of  the  N2  data  obtained.  An  asymmetry  of  this 

magnitude  will  not  have  an  appreciable  effect  on  the  llnewldth  determinations, 
but  could  produce  significant  systematic  errors  In  the  llneshlft 
determinations  because  the  llneshlfts  are  so  small.  This  asymmetry  Is 
discussed  later  In  Section  C. 

For  densities  larger  than  0.004  amagat  (~  3  torr  at  298  K) ,  the  rota¬ 
tional  Raman  lines  of  N2  ere  collision  broadened  giving  a  Lorentzlan  llneshape 
function.  Hence,  the  llnewldth  will  be  a  linear  function  of  density.  We  have 
therefore  taken  the  Lorentzlan  component  resulting  from  the  data  fitting 
procedure  described  In  the  previous  section  and  determined  the  density 
broadening  coefficients  B  for  each  temperature  studied  by  fitting  the 
llnewldths  to  a  straight  line.  The  data  and  the  linear  fit  (solid  line)  for 
the  S(10)  transition  are  given  In  Fig.  3  for  temperatures  of  298,  195,  and  80 
K.  Individual  points  Indicate  the  data,  and  the  solid  lines  represent  the 
least  squares  fits.  Uncertainties  (one  standard  deviation)  for  the  broadening 
coefficients  determined  from  the  fits  In  Fig.  3  are  1-2%. 

Table  I  summarizes  all  of  the  N2  self-broadening  coefficients  measured  In 
this  study.  These  values  are  graphically  Illustrated  In  Fig.  4  as  a  function 
of  rotational  quantum  number.  Odd  J  transitions  in  N2  were  not  studied.  For 
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S(4)  through  S(12),  the  specified  uncertainties  are  one  standard  deviation  for 
the  least  squares  fit.  For  S(0),  S(2),  S(14),  and  S(16),  data  were  obtained 
at  only  one  density,  and  the  uncertainties  shown  for  these  transitions  repre¬ 
sent  the  error  in  the  llnewldth  determination  of  a  single  line  profile. 

3  6  7 

Earlier  room  temperature  measurements  by  different  authors  ’•  have 
shown  considerable  disagreement  (30Z).  Table  III  compares  these  earlier 
results  with  our  present  results.  It  can  be  seen  that  the  three  earlier 
studies  show  significant  differences,  while  our  current  measurement  is  consis¬ 
tent  with  that  of  Ref.  7.  The  results  of  the  three  previous  studies  that  are 
shown  in  Table  III  were  obtained  at  densities  where  adjacent  lines  appreciably 
overlapped  with  one  another,  whereas  the  present  results  were  obtained  at 
densities  where  there  was  no  overlapping.  We  postulate  that  the  differences 
in  the  experimental  results  of  the  previous  studies  given  in  Table  III  are  due 
to  systematic  errors  arising  from  the  overlapping  of  adjacent  lines. 

Theoretical  calculations  of  the  collision  broadened  llnewldth  Raman  have 
been  performed  by  Grey  and  Van  Kranendonk, Srlvastava  and  Zaldl,^^  and 
Robert  and  Bonomy.  The  difference  in  these  three  approaches  is  in  the 
treatment  of  the  short  range  Intermolecular  interactions.  The  results  of 
these  different  approaches  are  compared  with  our  measurements  in  Fig.  3,  where 
the  room  temperature  broadening  coefficients  are  shown  as  a  function  of  rota¬ 
tional  quantum  number  J.  We  also  show  the  experimental  values  of  Janrniu  et 

3  11 

al.  for  comparison.  It  should  be  noted  that  Srlvastava  and  Zaldl  have  used 

O 

an  adjustable  parameter  to  obtain  their  fit  to  the  data  of  Jammu  et.  al., 
whereas  the  calculations  of  Grey  and  Van  Kranendonk^®  and  Robert  Bonomy^ “  do 
not  use  any  adjustable  parameters.  Fig.  5  shows  the  ab  initio  calculations  of 
Robert  and  Bonamy*  to  be  in  the  best  agreement  with  our  measurements. 

Previous  to  this  work,  only  atomic  foreign-gas  (He  and  Ar)  broadening  of 
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the  N2  transitions  had  been  Investigated.  We  have  also  measured  diatomic 
foreign-gas  (O2)  broadening  coefficients  at  195  K  and  295  K  for  the  four 
strongest  N2  transitions.  (Xir  O2  results  in  Table  II  show  that  the  foreign-gas 
broadening  by  ©2  Is  typically  85-90%  of  the  self-broadening  at  both  tempera¬ 
tures.  Thus  the  foreign-gas  broadening  by  0^  is  approximately  30%  larger  than 
the  foreign-gas  broadening  by  the  atomic  gases. 

B.  Line  Broadening  in  H2 

The  Raman  llnewldth  for  H2  has  a  much  more  complex  density  dependence 
than  that  for  N2.  This  results  mainly  from  the  fact  that  the  H2  density 
broadening  coefficient  is  approximately  a  factor  of  30  smaller  than  the  N2 
broadening  coefficient,  and  thus  Dlcke^^  or  collisional  narrowing  of  the  Raman 
llnewldth  can  be  observed^^* in  the  low  density  regime  if  the  resolution  is 
high  enough.  A  summary  of  the  relevant  collisional  processes  that  contribute 
to  the  llnewldth  at  various  densities  is  given  in  Ref.  14. 

For  densities  above  a  certain  cutoff  density  (p^  defined  below),  the 
Raman  llnshape  function  is  Lorentzlan  with  a  FWHM  llnewldth  dv  that  can  be 
expressed  as 


Av  =  A/p  +  Bp  (  1 ) 

where  A(MHz  amagat)  is  a  coefficient  proportional  to  the  self-diffusion 
coefficient  D^Ccm^  amagat  sec  ^).  For  forward  scattering, A  is  equal  to 
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where  is  the  Raman  transition  frequency  in  cm“^. 

This  model  is  known  as  the  diffusion  model  of  the  Raman  linewidth  and  it 
diverges  as  the  density  approaches  zero.  There  is  a  cutoff  density,  p^.,  at 
which  this  model  predicts  a  linewidth  that  is  10%  larger  than  that  predicted 
by  the  "hard  collision"  llneshape  theory. We  have  previously  determined 
that 

=  3.33  A/Avp  (3) 

where  AVjj(MHz)  is  the  FWHM  Doppler  width  of  the  Raman  transition  given  by 
(forward  scattering) 

AVjj  «  7.15x10"^  Vr  (T/m)^''^ 

where  is  in  MHz,  T  is  the  temperature  in  K  and  m  is  the  mass  in  amu.  For 
densities  larger  than  p^,  the  llneshape  is  Lorentzian  and  the  widths  can  be 
effectively  modelled  using  Eq.  1. 

For  example,  the  experimental  data  set  for  the  S(l)  transition  at  298  K 
is  given  in  Fig.  6.  From  the  Voigt  fit,  we  determine  both  the  Lorentzian  and 
Gaussian  contributions  to  the  experimental  linewidth.  Since  the  Raman 
llneshape  is  Lorentzian  for  all  densities  shown  in  Fig.  6,  the  Gaussian 
contribution  is  due  to  the  spectral  width  of  the  pulsed  dye  laser.  The 
Lorentzian  and  Gaussian  contributions  to  the  linewidth  are  plotted  in  Fig.  6, 
where  we  see  that  the  Gaussian  contribution  is  approximately  constant  (between 


80-100  HHz)  while  the  Lorentzlan  contribution  varies  approximately  linear  In 
density. 

Unfortunately,  the  100  MHz  resolution  of  our  current  laser  sytem  does  not 
allow  the  Dlcke  narrowing  of  the  H2  rotational  Raman  llnewldth  to  be  resolved 
with  any  precision.  We  have  therefore  chosen  a  fitting  procedure  that  uses 
determined  from  our  previous  study^^  to  calculate  A  from  Eq.  2.  The 
Lorentzlan  component  of  the  experimental  llnewldth  for  densities  larger  than 
Pj.  Is  then  fit  to  Eq.  1  with  one  free  parameter,  B.  The  parameters  used  In 
this  fitting  procedure  are  given  In  Table  TJ . 

The  density  broadening  coefficients  B  determined  using  this  procedure  are 
given  in  Table  V,  where  the  uncertainties  shown  are  one  standard  deviation  for 
the  one  parameter  linear  least  squares  fit.  This  fit  Is  Illustrated  by  the 
solid  line  In  Fig.  6  for  the  S(l)  transition.  The  maximum  correction  to  the 
llnewldth  due  to  Including  the  collislonal  narrowing  term  in  Eq.  1  Is  on  the 
order  of. 5-10%. 

We  have  also  Included  In  Table  V  data  from  the  previous  work  of  Van  Den 
Hout  et.  al.^^  and  Cooper,  May,  and  Gupta. There  are  small  discrepancies 
with  the  results  of  Refs.  16  and  18.  Some  of  these  differences  could  be  due 
to  small  errors  that  arise  In  our  estimate  of  A  for  these  rotational 
transitions.  We  emphasize  that  the  N2  results  of  Section  III.B  are  free  of 
this  problem  since  Dlcke  narrowing  Is  negligible. 

C.  Line  Shifts  In  N2  and  H2 

In  all  of  the  above  cases  for  which  density  broadening  measurements  were 
reported,  we  have  also  made  density  llneshlft  measurements.  In  all  cases  the 
llneshlfts  were  only  a  few  percent  of  the  llnewldth,  and  hence  any  asymmetries 
In  the  lineshapes  need  to  be  carefully  analyzed. 


The  asymmetry  In  FIr.  2  Is  typical  of  that  observed  in  our  data  for  the 
higher  density  N2  data.  The  lower  density  N2  data  contained  an  as3nnmetry  that 
occurred  only  further  out  in  the  wing  of  the  spectral  profiles.  About  802  of 
the  total  number  of  line  profiles  at  all  densities  and  temperatures  exhibited 
this  small  asymmetry,  where  t'^e  data  were  slightly  larger  than  the  symmetric 
fit  on  the  high  frequency  side  of  the  line.  Roughly  202  of  the  line  profiles 
were  symmetric  and  only  rarely  did  the  asymmetry  appear  such  that  the  data 
were  smaller  than  the  fit  on  the  high  frequency  side.  It  was  determined  that 
this  asymmetry  was  not  due  to  the  saturation  of  the  detection  system  by  either 
the  probe  beam  or  the  Raman  signal.  By  scanning  the  dye  laser  backwards,  it 
was  determined  that  the  asymmetry  always  remains  such  that  the  data  are  larger 
than  the  fit  on  the  high  frequency  side  of  the  line,  where  the  high  frequency 
side  corresponds  to  a  larger  Raman  shift.  For  densities  above  0.5  amagat  and 
for  the  pump  laser  intensities  used  here,  this  asymmetry  is  not  consistent 
with  the  measured  magnitudes  of  A.  C.  Stark  effects  in  N2.^^  At  this  time  we 
are  unable  to  determine  the  cause  of  these  asymmetries. 

In  N2,  both  the  self-density  shifts  and  the  O2  foreign-gas  density  shifts 
were  found  to  lie  in  the  range  50-150  MHz/amagat  for  the  transitions  and 
temperatures  discussed  above.  There  was  no  uniform  dependence  with  either 
temperature  or  quantum  number  J.  More  importantly,  there  were  inconsistencies 
observed  in  both  the  data  acquired  during  the  course  of  a  single  day  and  in 
the  data  collected  from  day-to-day.  These  inconsistencies  have  forced  us  to 
conclude  that  the  line  shift  measurements  in  N2  are  substantially  disturbed  by 
the  asymmetric  lineshapes.  The  data  can,  however,  be  used  to  establish  an 
upper  limit  to  the  magnitude  of  the  line  shifts.  For  all  transitions  and 
temperatures  Indicated  for  N2  in  Tables  I  and  II,  the  density  llneshifts  are 
less  than  150  MHz/amagat.  Future  elimination  of  the  asymmetry  would  enable  us 


to  detect  llneshlfts  as  small  as  10  MHz/amagat  with  our  current  experimental 
apparatus. 

The  density  shifts  observed  In  H2»  unlike  those  In  N2,  were  free  of  the 
Inconsistencies  mentioned  above.  Data  for  S(l)  In  H2  at  295  K  are  Illustrated 
In  Fig.  7,  where  the  density  shift  Is  defined  as 

fiv  ■  ^5) 

where  vr(p)  Is  the  Raman  transition  frequency  at  the  density  p.  Hence  a  nega¬ 
tive  shift  results  In  a  smaller  Raman  frequency.  Density  shifts  were  plotted 
as  a  function  of  density  difference  between  the  sample  and  reference  cells 
rather  than  density  because  data  were  collected  for  reference  cell  pressures 
of  50  torr,  100  torr,  and  200  torr.  The  temperature  Is  295  K  for  the  data  of 
Fig.  7.  A  summary  of  our  H2  density  shift  results  is  given  in  Table  VI,  along 
with  the  results  from  Refs.  16  and  18.  The  agreement  between  our  results  and 
those  of  Ref.  18  Is  consistent  with  the  combined  uncertainties  for  the  two 
experiments. 

It  should  be  noted  that  the  density  shift  for  the  rotational  Raman  line 
In  H2  Is  substantially  different  from  that  observed  for  the  vibrational  Raman 
lines. For  example,  at  80K,  the  shift  of  the  0(1)  vibrational  transition  is 
13  times  larger  than  the  S(l)  rotational  Raman  line.  This  difference  results 
from  the  perturbation  of  the  vibrational  frequency  of  the  H2  molecule  during 
an  elastic  collision;  an  effect  not  present  for  the  rotational  transition. 


IV.  DISCUSSION 


A  simple  qualitative  understanding  of  this  temperature  dependence  for  the 

N2  rotational  Raman  broadening  data  given  In  Table  I  can  be  obtained  from  a 
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consideration  of  Anderson's  Impact  theory  of  the  colllslonal  broadening. 

The  linewidth  can  be  expressed  as 

Av  *  n  <va(v)>  =  n  <v><a(v)>  (6) 

where  n  is  the  molecular  number  density  (cm~^),  v  is  the  relative  velocity  of 
the  collision  partners,  a  Is  the  colllslonal  broadening  cross  section,  and  the 
brackets  denote  an  average  over  the  velocity  distribution  of  the  gas.  As  can 
be  seen  from  Eq.  6,  the  temperature  dependence  of  the  density  broadening 
coefficient  is  Implicit  in  the  velocity  dependence  of  the  colllslonal  cross 
section.  For  example,  if  we  assume  a  "hard  sphere”  crodel  for  the  cross 
section  (velocity  Independent),  the  density  broadening  coefficient  (Av/n) 
measured  in  our  experiment  would  be  proportional  to  T®*^. 

The  broadening  cross  section  is  composed  of  contributions  from  several 
different  type  of  collisions.  Both  elastic  phase  perturbing  collisions,  and 
rotationally  inelastic  collisions  contribute  to  the  linewidth.  However,  there 
are  several  cases  where  the  contribution  from  elastic  collisions  is  small. 

For  example,  the  llnewidths  for  isotropic  vibrational  0  branch  Raman  scatter¬ 
ing  have  no  contribution  from  elastic  collisions  if  collisions  perturb  the 
ground  and  final  vibrational  level  the  same  amount.  This  appears  to  be  a  good 
approximation  for  all  cases  except  hydrogen^^. 

Rotational  Raman  llnewidths,  on  the  other  hand,  have  contributions  from 
phase  perturbing  elastic  collisions.  This  is  the  primary  reason  why  rotational 


Raman  transitions  usually  have  larger  llnewldths  than  vibrational  Raman 
transitions.  However,  for  the  case  of  N2,  CO,  and  CO2,  it  has  been 
calculated^*^  that  elastic  collisions  contribute  no  more  than  15%  of  the  cross 
section.  Therefore,  in  most  cases,  the  primary  problem  is  to  calculate  the 
velocity  dependence  of  the  cross  section  for  rotatlonally  inelastic 
collisions. 

The  calculation  of  the  colllslonal  cross  section  has  been  the  subject  of 
considerable  research.  A  good  review  of  the  various  approaches  to  this  calcu¬ 
lation  for  Raman  transitions  has  been  given  by  Srivastava  and  Zaldl^^.  One 

approach  to  this  calculation  is  to  expand  the  long  range  part  of  the  inter- 
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molecular  potential  in  a  multipole  expansion  and  then  calculate  the 

contribution  to  the  cross  section  of  each  term  in  this  expansion.  Using  this 
approach.  Gray  and  Van  Kranendonk^®  have  shown  that  for  molecules  such  as  N2, 
CO,  and  CO2,  the  most  important  terms  in  this  expansion  for  self-broadening 
are  the  quadrapole-quadrapole  (OQ)  potential  and  the  dispersion  potential. 

IQ 

It  can  be  shown  ’  in  general  that  the  velocity  dependence  of  these 
contributions  to  the  cross  section  scale  as 

<a>  «  <v~2/(n-l)> 

where  n  is  the  exponent  of  the  radial  dependence  of  the  long  range  part  of  the 
potential  in  the  form  V(r)  ~  r~”.  For  the  QQ  potential,  n  *  5,  while  for  the 
dispersion  potential  n  ■  6.  Subsituting  Eq.  7  into  Eq.  6  and  remembering  that 
V  a  T*^,  we  find  that  the  density  broadening  coefficient  should  scale  with 


temperature  as 


(Av/n)  «  T 


.(n-3)/2(n-l) 


(8) 


I 
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From  this  simple  argument,  we  would  expect  our  density  broadening  coefficients 
to  scale  as  i®*25-0.3^  depending  on  the  relative  contribution  to  the  cross 
section  of  the  QQ  and  dispersion  forces. 

To  test  this  simple  model  of  the  temperature  dependence,  we  have  fit  our 
data  In  Table  I  to  the  equation 


B  -  Bo(T/295)^  (9) 

where  B^  Is  the  density  broadening  coefficient  at  295  K.  The  results  are 
given  In  Table  VII  as  a  function  of  Initial  rotational  state.  We  find  that 
0.26  <  Y  <  0.39,  In  reasonable  agreement  with  this  simple  picture.  We  believe 
that  most  of  the  scatter  In  y  as  a  function  of  J  Is  due  to  the  fact  that  a 
much  larger  temperature  range  needs  to  be  Investigated  to  accurately  determine 
this  small  temperature  dependence* 

However,  the  scaling  given  In  Eq.  9  may  not  be  valid  over  a  wide  tem¬ 
perature  range.  By  theoretically  examining  temperature  dependence  of  the 
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rotational  Raman  llnewldth  over  a  wide  range.  Pack  proposes  a  quadratic 
temperature  scaling  law  for  density  broadening  coefficients  that  can  be 
written  as 


B  -  Bq  +  Cj,(T-295)  +  Do(T-295)2 


(10) 


where  6^  Is  the  broadening  coefficient  at  298  K  (same  as  In  Eq.  8),  and  and 
Dq  are  fitting  coefficients.  We  have  previously  used  this  temperature  scaling 
formula  In  modeling  Raman  llnewidths  in  H2.^^  The  data  from  Table  I  has  been 
fit  to  Eq.  10  assuming  Dq  »  0,  and  the  results  are  also  given  In  Table  VII. 

We  find  that  this  expression  Is  just  as  useful  as  Eq.  9  in  modelling  the 
dependence  of  the  density  broadening  coefficient  over  the  limited  temperature 
range  Investigated  In  our  experiment.  Both  Eqs.  9  and  10  can  be  used  to 
predict  B  to  better  than  3Z  over  the  range  of  80-300  K.  This  result 
emphasizes  the  need  for  further  measurements  at  high  temperature  to  accurately 
determine  whether  an  exponential  or  quadratic  scaling  law  Is  approplate  to 
model  the  Raman  llnewidths  over  a  large  temperature  range. 

There  are  few  molecular  systems  for  which  the  rotational  Raman  llnewidths 
have  been  determined  as  a  function  of  temperature,  and  no  other  studies  in 
N2«  However,  there  have  been  extensive  studies  of  the  temperature  dependence 
of  absorption  llnewidths  in  and  C02.^^  The  absolute  magnitude  of  the 

llnewidths  determined  from  absorption  will  be  different  from  those  determined 
from  Raman  scattering.  This  Is  due  to  the  fact  that  different  states  are 
involved  in  the  calculation  of  the  collision  cross  section.  However,  the 
temperature  dependence  will  be  similar  If  the  same  terms  In  the  multipole 
potential  are  responsible  for  the  broadening.  This  is  the  case^®  for  both  CO 
and  CO2  where  the  most  Important  terms  in  the  potential  are  the  same  as  for 
N2« 

Most  of  the  data  for  llnewidths  in  the  literature  is  given  in  terms  of 
pressure  broadening  coefficients.  We  can  compare  results  by  noting  that  y  in 
Eq.  9  is  related  to  n  by  y  ■  1  -  n,  where  n  Is  the  exponent  of  the  temperature 
for  pressure  broadening  coefficients.  We  find  that  y  for  CO  and  CO2 
absorption  llnewidths  ranges  cover  a  similar  range  of  values  as  observed  for 
the  rotational  Raman  llnewidths  for  N2. 


The  details  of  the  exact  theoretical  temperature  dependence  of  the  rota¬ 


tional  Raman  llnewldths  requires  an  ab  Initio  calculation  of  the  velocity 
dependence  colllslonal  cross  section,  and  Is  beyond  the  scope  of  this  study. 
However,  the  precision  of  the  data  presented  here  should  allow  a  detailed 
comparison  with  llnewldths  calculated  from  the  various  theoretical  approaches 
thus  allowing  new  theoretical  advances  to  be  made. 


V.  CONCLUSIONS 


The  density  self-broadening  coefficients  have  been  measured  at  tempera¬ 
tures  of  80  K,  195  K,  and  295  R  for  the  even-J  rotational  Raman  transitions  in 
N2.  The  temperature  dependence  of  these  coefficients  was  found  to  fit  both 
linear  and  T^  models  equally  well.  The  room  temperature  results  of  this  work 
are  in  agreement  (less  than  5%  differences  for  J  >  2)  with  the  most  recent  ab 
Initio  calculations.  Foreign-gas  broadening  of  these  N2  transitions,  due  to 
O2,  was  also  measured  and  found  to  be  10-15Z  smaller  than  the  self-broadening 
coefficients.  We  have  also  determined  that  the  density  self-shifts  and  O2 
foreign-gas  shifts  for  these  same  N2  lines  are  no  larger  than  150 
MHz/amagat.  In  a  previous  letter^,  we  have  used  the  above  results  In  the 
modeling  of  the  temperature  dependence  of  the  steady-state  Raman  gain 
coefficient  In  N2. 

The  rotational  Raman  density  broadening  and  shifts  In  H2  were  measured  at 
80  K  and  295  K.  Our  results  were  In  general  agreement  with  previous  experi¬ 
mental  work.  This  agreement  of  our  H2  data  with  other  independent  measure¬ 
ments  provides  confidence  In  the  new  temperature-dependent  N2  results  pre¬ 
sented  above. 
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FIGURE  CAPTIONS 


Figure  1.  Schematic  of  experimental  apparatus  for  quasi-cw  stimulated  Raman 
spectroscopy* 

Figure  2.  Example  of  S(10)  llneshape  in  N2  at  195  K  for  pressures  of  10  torr 
(narrower  profile)  and  400  torr.  Individual  dots  are  data  and  the 
solid  lines  are  Voigt  fits. 

Figure  3.  S(10)  llnewidths  in  N2  at  temperatures  of  80  K  (□),195  K(0)» 

and  295  K  (A).  Solid  lines  are  linear  least  square  fits  to  the 
data. 

Figure  4.  J  dependence  of  the  N2  rotational  Raman  broadening  coefficients  for 
80  K  (□),  195  K  (O).  and  295  K(A). 

Figure  5.  Comparison  of  theoretical  calculations  with  experimental  data  for 
colllslonal  broadening  of  the  rotational  Raman  lines  in  N2  at  room 
temperature.  The  calculations  are  from  Ref.  10  (long  dash).  Ref. 

11  (short  dash),  and  Ref.  12  (solid  line).  The  experimental  data 
are  from  Ref.  3  (□)  and  this  work  (A). 

Figure  6.  Lorentzian  (  □  )  and  Gaussian  (A)  components  from  Voigt  fits  to 
data  as  a  function  of  density  for  S(l)  in  82*  Solid  line  is  a 
linear  one  parameter  fit  to  Eq.  I. 

Figure  7.  Example  of  data  for  llneshift  determination  for  S(l)  in  H2  at 
295  K.  The  solid  line  is  the  least  squares  fit. 


B  (MHz/amagat) 


295  K 


195  K 


S(J) 


295  K 


195  K 


6  2970  ±  130  2690  ±  120 
8  2770  ±  80  2600  ±  100 
10  2750  ±  50  2410  ±  100 
12  2450  ±  70  1980  ±  100 
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TABLE  IV.  H2  rotational  Raman  parameters  for  an  ortho:para  ratio  of  3:1 


TABLE  VII.  Fit  parameters  for  temperature  dependence  of  density 
self-broadening  coefficients  in  N2* 
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